Mean-field dynamics of a quantum dot - microcavity system 
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Mean-field evolution equations for the exciton and photon populations and polarizations (Bloch- 
Lamb equations) are written and numerically solved in order to describe the dynamics of electronic 
states in a quantum dot coupled to the photon field of a microcavity. The equations account for 
phase space filling effects and Coulomb interactions among carriers, and include also (in a phe- 
nomenological way) incoherent pumping of the quantum dot, photon losses through the microcavity 
mirrors, and electron-hole population decay due to spontaneous emission of the dot. When the dot 
may support more than one electron-hole pair, asymptotic oscillatory states, with periods between 
0.5 and 1.5 ps, are found almost for any values of the system parameters. 
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I. INTRODUCTION 

The coupling of light with an atomic state depends 
basically on two factors: (a) the spread of the wave func- 
tion or dipole moment of the atomic state, and (b) the 
density of photon states in the energy region of interest. 
In the search for optimizing this coupling, atomic Ry- 
dberg states and three-dimensional microcavities, where 
the photon modes exhibit a discrete spectrum, have been 
used. One of the sahent realizations is the construction 
of micromasersi. 

More recently, excitonic states in quantum dots have 
proven to be excellent candidates to replace atomic Ry- 
dberg states in a microcavity. From one side, quantum 
dots can easily be incorporated inside the microcavity by 
means of modern semiconductor technology, thus avoid- 
ing the complicated step in which a beam of atoms in Ry- 
dberg states is prepared. On the other hand, the effective 
Bohr radius in a semiconductor is a factor of woe/ (meo) 
larger than the corresponding radius in vacuum, where 
mo and eq are, respectively, the electron mass and dielec- 
tric constant of vacuum, and m, e are the corresponding 
magnitudes in the semiconductor. In GaAs, for example, 
this factor is around 200. As a result, the ground-state 
exciton wave function is as spread as an atomic Rydberg 
state. Transitions energies are typically in the infrared 
range. 

In the present paper, we study the dynamics of a quan- 
tum dot-microcavity system, focusing on the stationary 
or asymptotic photonic states in the cavity. A previ- 
ous, very serious, study of this dynamics^ has been mo- 
tivated by recent experimental work'^'^'^, and theoretical 
papersSii, in which the exciton-cavity coupling is charac- 
terized, and the system is shown to serve as an efficient 
source of single triggered photons, and a possible source 
of entangled photons. 

The results of Ref. |2 rely on the assumption that the 
quantum dot may support only one excitonic state (a 
shallow dot). We relax this assumption, allowing more 
than one excitonic state in the dot. Coulomb interac- 
tions among particles are shown to play an important 



role in the dynamics, shifting the single-pair levels, mod- 
ifying the Rabi frequencies, and causing interference ef- 
fects. The most relevant consequence of including addi- 
tional pair states in the dot, however, is shown to be the 
existence of asymptotic oscillatory states, in which the 
number of photons oscillates around a mean value, with 
periods between 0.5 and 1.5 ps (for the set of parameters 
used in the calculations). We keep the assumption that 
the excitonic states are coupled to a single photon mode. 
This means that the separation between photon modes 
in the cavity should be larger than the splitting between 
excitonic states in the dot. As the latter magnitude is 
of the order of a few meV, our results will be valid for 
cavities with a radius of 0.5 /j,m or smaller—. 

A second, simplifying, working hypothesis is the mean 
field approximation. For electrons, it means the neglect 
of correlation terms in the evolution (Bloch) equations. 
For light, we use a semiclassical or coherent (Lamb) de- 
scription, in which the field is characterized by an am- 
plitude and a phase. Let us stress that, in the intervals 
of parameters we consider, when the number of photons 
proves to be significantly different from zero the statis- 
tics of the photon field is Poissonian or near Poissonian, 
at least in the onc-exciton dot^, in such a way that the 
coherent approximation for the photon field is justified. 
The reward from this simplification is that the system of 
evolution equations is relatively small, allowing a qualita- 
tive analysis of its solutions, and the extension to deeper 
dots, which may support more than one electron-hole 
states. 

The plan of the paper is as follows. In Sec. m the basic 
evolution (Bloch-Lamb) equations are obtained, and the 
initial conditions used to solve them are discussed. Sec. 
mil is devoted to the presentation of results. The one-, 
two-, and three-states dots coupled to a single photon 
mode (with a given polarization) are studied. Addition- 
ally, the two-states dot coupled to a polarization degen- 
erated light mode is also considered. Finally, concluding 
remarks are given in Sec. IIVI 
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II. THE THEORY 
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Our model microcavity is assumed to have an ideal 
cylindrical shape. Its radius is lower than 0.5 /im, in 
such a way that the spacing between photon modes is 
of the order of 20 meVi^ We will be interested in one 
of these modes. To fix ideas, we speak about the lowest 
mode, polarization degenerated, which is called HEn ac- 
cording to the terminology for modes in a guided. In the 
paper, we will conventionally denote the two polarization 
components as cr(+) and cr(_). 

A quantum dot is located at a point along the cavity 
axis in which the electric field has a maximum. We as- 
sume a GaAs dot. The dot has also a cylindrical shape. 
Its height is lower than 10 nm, in such a way that light- 
hole sub-bands have higher energy than the heavy-hole 
ones and, in a first approximation, could be ignored in the 
dynamics. A second reason for not including light holes 
is that their coupling to the photon modes is a factor of 
3 less intense than the heavy hole couplingiii. 

Electron-hole pairing is induced by both Coulomb in- 
teractions and the coupling to the photon mode. Thus, 
in addition to the electron, hole and photon populations: 
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In these equations, n and n label electron and hole 
single-particle states in the dot. ejj and h\ are the corre- 
sponding creation operators. a(+) and a(_) are operators 
for the two photon polarizations. In the k functions, we 
have explicitly indicated the electron and hole states that 
are coupled to a polarized mode. For example, in k'^"''' 
the coupling of a spin-down electron state (rrij — —1/2) 
with a "spin-up" hole state (coming from a rrij = —3/2 
electron state in the valence band) is due to the CT(+) 
mode, and it is reinforced by Coulomb interactions. 

Notice also that we have assumed pairing between spe- 
cific electron and hole states, n and n (orbitals). We will 
use a harmonic oscillator basis. Thus, for a given n char- 
acterized by radial and angular momentum projection 
quantum numbers, n = (fc, I), the hole state coupled to n 
is fi = (fc, —I). The pair is created in a zero total angular 
momentum state, corresponding to the selection rule for 
interband transitionsiS. 

Mean values, (...), come from averaging with a density 
matrix. We are interested in the very low temperature, 
T — > limit. This density matrix is time dependent. 
Thus, the population and polarization functions, Eqs. 
(|1I6|I . are time-dependent functions. We will derive dy- 
namical equations for their time evolution (Bloch-Lamb 
equations). 

A further simplification equates pnn — p'nl — Pn- This 
means that the dot is initially neutral, and that we are 
interested in the dynamics over a time interval smaller 
than (or of the order of) the relaxation times in the dot 
(the coherent regimaii). In fact, we will follow the dy- 
namics in the time interval (0, 20 ps). 

The Hamiltonian describing the quantum dot - micro- 
cavity system can be written in the following form: 
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The En functions refer to the single-particle energies of 
the used basis states. On the other hand, tkn are the ma- 
trix elements of the external (quantum dot) confinement 
potential. (3 is the strenght of Coulomb interactions, and 
g the strenght of the electron-photon coupling (assumed 
constant). The energy of the photon mode is written as 



A. Bloch-Lamb equations 

The dynamical equations are obtained by taking time 
derivatives of the occupation and polarization functions 
pi6|l . One gets for p„, for example: 



ih^ = {[eie^,H]). (8) 
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Once the commutator is computed, we take the mean- 
field or quasiparticle contribution (in which mean values 
of products of four or more operators are expressed in 
terms of products of occupations and polarization func- 
tions), and neglect the collision termsii. The result for 
Pni is: 



dt 



^ Im I K~^i*^{n,j\l/r\j,n)K-. 
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and a similar equation for p„|. Notice that Coulomb 
interactions preserve spin. Thus, if n is a spin-down state 
(as assumed in Eq. JHJl), the j states entering the sum 
should be spin-down states also. This is the reason why 
only functions enter the sum. 

The equation for the photon number is straightfor- 
wardly obtained also, and leads to: 
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Notice the conservation of the polariton number, as 
it follows from Eqs. © and H1U|) (or directly from the 

Hamihonian, Eq. 0): dN^^l^} /dt = 0, where: 
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The equations for the polarization functions, k, are ob- 
tained in the same way. We shall stress, however, that 
the modulus, is determined from Eq. 0, and it fol- 
lows that \K-nf^\ — ^ Pni(l — Pni)- Nevertheless, we are 
forced to use the equation for k'^^ in order to determine 
its phase. 

In the same way, the modulus of (T(_|_-) (understood as 
(a(+))) is determined by Eq. IjlUI) . but its phase is not. 
Thus, we must write the equation for cr(^). 

Both a and k should contain a rapidly varying phase 
factor, e~*"*, due simply to the frequency of rotation of 
the electric field. For example, w ^ 2 x 10'^ ps^^ for the 
GaAs band gap. In contrast, we expect characteristic 
frequencies of the order oi g/h ^ 1 ps~^ in the variation 
of pn or |(t|. Thus, we will use a kind of Rotating Wave 
Approximation, and write explicitly the trivial (rapid) 
phase dependence of a and k: 
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etc. By definition, S(+) — |cr(+)|, and n'it^ 
equations for cj) and (pn are, then: 
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and similar equations for cf)^^^ and (/)„ . We may easily 
identify in Eqs. (|14I15|) the space-filling factors, 2p„ — 1, 
mean-field contributions to the pair energies, detuning 
factors, etc. 



Eqs. H9I15|I are Bloch equations for the electronic oc- 
cupations and polarizations in the quantum dot^^, in a 
mean field approximation and under a coherent regime, 
whereas (|1UI14II are the analogues of the semiclassical 
Lamb equations for a laser^^. Thus, we will call our 
system of evolution equations Bloch-Lamb equations. 
Losses and incoherent pumping are introduced in the 
next subsection. 



B. Phenomenological account of pumping and 
losses 

The main processes of the interaction between the 
quantum dot - microcavity system and the external 
medium are described in Ref. i3, and modeled by means 
of additional terms in the Liouville equation for the den- 
sity matrix. 

One of such processes is the decay of cavity photons 
by escaping through the cavity mirrors. We may account 
for this process by introducing a decay term, — >r|CT(+)p, 
into Eq. (|10|l . with a phenomenological decay rate, x. 
Roughly, we have h ~ hw /Q, where the quality factor of 
the cavity is in the range between 1000 and 5000 J> 
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In terms of the variables introduced in Eqs. (|12I13|I . 
the modified Eq. (|10|) is, thus, rewritten: 

^=fE4"^-(^^"^-^r')--wA (16) 

On the other hand, the electronic population, 



is modified by a continuous incoherent pumping of the 
quantum dot at a rate P(+), and the decay of electron- 
hole pairs by spontaneous emission into leaky modes at 
a rate 7. In Eq. ©, one should add the following terms: 
P(+)(l — Pni) ~ "fPni, where Pauli blocking is explicitly 
introduced in the pumping contribution. The resulting 
equation for pni is, then: 



dt 
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The set of equations 114I17II determine the magnitudes 
Pni, S(+), (?!)(+\ and n is obtained in terms of 

Pni as Pnj,(l — Pni)- Noticc that the phase of is 
not affected by the decay processes in this approxima- 
tion. The modulus of decays at a rate which is 
determined by the decay of pni- This is the analogue of 
the relation T2 = 2Ti between the decay times of pop- 
ulations and coherences, valid both in Atomic Physics^- 
and, in a certain way, in bulk semiconductors^^. Let us 
stress also that the r.h.s. of Eqs. H14I17|I depend only on 
differences of angles. Thus, we may eliminate one of the 
angles, for example (j), by introducing the new variables 

9n = 4>- 4>n- 

A formally similar set of equations may be written for 

S(-), ip^^K and (pn Notice that, in our approxima- 
tion, where Coulomb collision terms and other relaxation 
processes (due to the interaction with phonons, impu- 
rities, etc) are neglected, the relations pm = p^^n are 
preserved at all times, and p„i evolves independently of 

Pn]- 
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where N states is the number of single-particle states avail- 
able in the dot, and 
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n and j in Eqs. (|23I24|) represent spin-down states. Eqs. 
(|18l24f) are to be used in the numerical integration of 

Eqs. rmm . 



C. The t —> asymptotics 



Eqs. H14I17|I are to be solved with null initial condi- 
tions, i.e., there are no pairs in the quantum dot and no 
photons in the cavity at t 0. 

According to Eq. H17I) . the populations p„| rise as 
P{+)t and, consequently, k^^-* ~ y^P(^^)t. However, the 
small-^ asymptotics of the magnitudes S(4_), 0*^+-' and 
(f)^^'^ should be analysed with care because of terms like 
Kri^V'*(+): 6tc entering the equations. We found that the 
consistent behavior in the t ^ limit is the following: 



III. NUMERICAL RESULTS 

In order to perform the numerical integration of Eqs. 
(|14ll7f) . we shall give a precise meaning to the magni- 
tudes entering them. We are interested in a quantum dot 
which supports a finite number of single-particle states, 
N states- Thus, N states wiU bc fixcd by hand. A magnetic 
field, B = 7 Teslas, acts on the charge carriers in the dot. 
Both, the dot confinement and the magnetic field may be 
thought of as control parameters. Energy magnitudes are 
converted to angular frequencies by dividing by h. 

We will use parameters for GaAs. The single-particle 
energies are defined as: 



Pni = P{ + )t + 
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where the cyclotron frequencies are uJc = 2.625 B ps ^, 
ivi''^ = 1.173 B ps"^ (and B in Teslas). 5"^ = ±1/2 are 
the spin projections along B. ge^B/{2h) — 4.397 x 10"'^ 
ps-\ and ghfiB/{2h) = -4.397 x 10"^ pg-i^ xhc + sign 
in Eq. (^5)1 corresponds to the electron. The fact that, 
in our simple model, energy levels come from a parabolic 
confinement and a magnetic field has not a decisive im- 
portance. What really matters is the position of the bare 
levels, and the effects of Coulomb interaction among car- 
riers, as will become evident in the next sections. 

On the other hand, the matrix elements of the 
(parabolic) confinement potential are written as: 

fe, ,2 

= ^(2fc„ + |;„| + l), (26) 

LUc 

where the dot frequency is loq = 4.558 ps^^. The 
electron-photon coupling, g, which determines the Rabi 
frequency, is fixed to g/h — 1.5 ps~^, and the strength of 
Coulomb interactions is P/h — 3.871\/B ps~^. 

Notice that we have not included the nominal gap en- 
ergy into the electron energy, Eq. (|25|l . A shifted mode 
frequency is defined according to A = w — Egap/h. 

Concerning the loss and pumping parameters intro- 
duced in Sec. Ill Bl we will set the rate of spontaneous 
decay, 7, to 0.1 ps^^. The equations are integrated for 
differents sets of P, x, and A. A variation of A should 
be understood as varying the cavity dimensions. In fact, 
B, loq and P are the only parameters which can be ex- 
ternally controlled for a given dot-microcavity system. 

Once we have precisely defined our equations, we may 
turn to the analysis of the simplest system, in which the 
dot may support only one pair. 
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FIG. 1: (Color online) Number of photons, dot occupation 
and polarization in the one-pair dot for P = 1 ps~^, x = 
0.5 ps~^. The photon energy is A = 4.5 ps~^ (maximum 
resonance) . 



A. A shallow dot with Nstates ~ 1 

This case was extensively studied in Ref. 0. The 
authors obtain and solve the equations for the density 
matrix in a Pock basis composed by the electron-hole 
vacuum, |G), the exciton ground state, \X), and the n- 
photon states, \n). They write equations for the time 
evolution of the populations, pon.Gn, PXn.Xm and the 
coherences, pcn.Xn-i- To solve them, they should trun- 
cate the photon Hilbert space. For n < 100, for example, 
they have 300 equations. This is an exact treatment, in 
which the authors may compute correlation functions for 
the photon field, the spectrum of the emitted light, etc. 

In our model, however, the Nstates = 1 case is de- 
scribed by only three variables, pi, s, and the phase dif- 
ference (j) — 4>i. The photon field is assumed coherent. 
In spite of these simplifications, most of the results pre- 
sented in paper ^] concerning the system dynamics are 
qualitatively and quantitatively reproduced by our equa- 
tions. In the present section, we focus on new results, 
uncovered in paper 0. 

The first feature we observe in our calculations is the 
presence of the three regimes, mentioned in under 



which the quantum dot - microcavity system may op- 
erate. They can be roughly classified according to the 
relation between pumping and losses: (i) P 3> >f, (ii) 
P and (iii) x:^ P. 

We show in Fig. ^ an example of calculations corre- 
sponding to the intermediate regime (ii), in which P = 1 
ps~^ and X = Q.b Y>s~^ ■ This example reveals oscillations 
in the transient interval until the asymptotic values are 
reached. The oscillations of the number of photons in the 
cavity are due to non compensation between the rates of 
creation of photons and losses, and even to photon ab- 
sorption by the dot (when sin(0 — ipi) < 0). 

The asymptotic number of photons, computed at tf — 
100 ps as a function of A for a high quality cavity {x = 
0.1) is shown in Fig. [3 The pair frequency (pair energy 
divided by h) in the present case is 

impair — yJ^i + -C^J + 111 + tfi J 

- ^(l,l|l/r|l,l)=4.5ps-i. (27) 

Extreme resonance corresponds to A = uipair- The de- 
pendence Nph vs A in Fig. {P = 7 ps~^) resembles 
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FIG. 2: The number of photons as a function of the photon 
energy for (a) P = 7 ps~^, >c — 0.1 ps~^, and (b) P = 1 ps~^, 
X = 0.1 ps~\ 
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FIG. 3: (Color online) The same as in Fig. [T] but for P = 1 

ps~^, yc — 0.1 ps^^, and A = 8 ps~^. 



a step function with a mean value of 30 photons inside 
the resonance interval, and zero outside. If the cavity is 
such that the actual A lies in the abruptly varying re- 
gion, small changes produced by variations in B or in wo 
could switch between the two "states" of the cavity. Al- 
though the step-like variation could be an artifact of our 
simplified equations, an abrupt variation is expected. 

Similar step-like variations are shown in Fig. I^j, which 
corresponds to P = 1 ps~^. In addition, outside the reso- 
nance interval there are asymptotic oscillating solutions, 
whose characteristic behavior is illustrated in Fig. Ofor 
A = 8 ps~^. Notice that, in the oscillating solution, the 
quantum dot is periodically emitting (sin((/i — 0i) > 0) 
and absorbing (sin(0— (^\) < 0) photons from the cavity 
with a period around 2.5 ps. The mean number of pho- 
tons is only around 0.1, thus our coherent approximation 
could be very rough. Such periodic solutions are, how- 
ever, very common in the larger systems with N states > 1, 
as will be seen below. 

From the point of view of the theory of nonlinear evo- 
lution differential equations^^, in the resonance interval 
the solution of the system H14I17|I is attracted by a stable 
critical point (an spiral) , whereas the asymptotic periodic 
solution outside the resonace interval for A corresponds 



to a limit cycle. This view reveals the existence of a sec- 
ond stable attractor (also an spiral) in the resonance in- 
terval, which is not seen in calculations because the used 
initial conditions is not in its basin of atraction. As we 
move out of the resonance interval, one of the stable crit- 
ical points dissapears (a bifurcation) and the limit cycle 
emerges. In addtition, at any value of A there is a third 
critical point of hyperbolic character (with one unstable 
direction). The existence of two asymptotically stable 
states of the quantum dot-microcavity system, how to 
reach them, and even how to reach the hyperbolic point 
by using control of chaos^^, the nature of the limit cycles, 
etc are very interesting questions which deserve further 
investigation. In the larger, N states > 1, systems the 
dynamics is expected to be even more complex. In the 
present paper, we restrict the analysis to the situation 
with null initial conditions, described in Sec. Ill CI and 
leave the more general questions for a later work. 



B. The Nstates ~ 2 dot coupled to CT(+) light 

In view of the fact that, in our equations, the CT(+)- 
polarized mode evolves independently of the tT(_) mode, 
the next nontrivial case corresponds to a dot which may 
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FIG. 4: The single-pair levels as a function of the photon 
energy for P = 1 ps~^, >r = 0.1 ps~^. 



support up to two pairs, Nstates = 2, interacting with 
the same photon mode. To be specific, we will consider 
coupling to the cr(+) mode. 

The dynamics, in many ways, has features very similar 
to those of the simpler one-pair system. In particular, the 
three operation regimes are observed. In addition, there 
are new characteristics, which are briefly outlined below. 

The first new (expected) feature is the two-peak struc- 
ture in the number of photons as a function of A as a 
result of individual resonances with the single-pair lev- 
els. Coulomb interactions are crucial in determining the 
position of these levels. Indeed, we have: 



' E + -J- E^^^ + -I- t^''^ 



n, n\l/r\n, n) 



2(3 



(28) 



In our model quantum dot, the pair frequencies at 
(3 = are igual to 17.32 and 20.99 ps~^. The account of 
Coulomb interactions move them to the interval from -2 
to 10 ps~^, as shown in Fig. ^ and make the single-pair 
levels dependent on the occupations pn- But Coulomb 
interactions modify also the phases and occupations (see 
Eqs. ifT^ and lfT7|l ) causing strong changes in the sta- 
tionary solutions. We illustrate in Fig. |S1 these effects 
for the system with parameters P = 1 ps~^ and k = 0.1 




FIG. 5: The number of photons as a function of the photon 
energy for the Nstates = 2 system with P = 1 ps~^, x = 
0.1 ps~^. In the upper panel. Coulomb interactions are not 
included (/? = 0). 



ps . 

A second interesting characteristic is that asymptotic 
oscillatory solutions are very common, in particular in- 
side the resonance intervals. It means that, as a rule, 
the number of photons and the occupations will vary pe- 
riodically around a mean value. In our Fig. |S1 where 
the number of photons is computed at fixed tf = 100 
ps, oscillatory solutions are seen as apparently random 
behavior of Nph- 

It is interesting also to look at the relative phases, (l>—(f>i 
in each of the resonance intervals. Let us consider, for 
example, the system described in Fig. ]Ejp. In the up- 
per interval, 3 < A < 10 ps^^, the two single-pair levels 
emit or absorb almost in phase (see Fig. Et)- In the 
lowest interval, — 5 < A < ps~^, one level is emitting, 
and the second one emits and absorb, but they do so "in 
phase" . That is, maximum emission occurs at the same 
time in both levels, and maximum absorption in one of 
them occurs at coincidence with minimum emission in 
the other. The situation is depicted in Fig. ^p. Finally, 
in the central peak around A « 2.5 ps^^, which is the re- 
sult of interference and dissapears when parameters are 
changed, emission and absorption in the two levels are 




FIG. 6; (Color online) sin(<^ — (pi) in the Nstates = 2 dot and 
cavity parameters P = 1 ps~^, m: = 0.1 ps~^, and (a) A = 7 
ps-\ (b) A = -3 ps'\ and (c) A = 2.5 ps'^ 



"in counterphase" , as can be seen in Fig. Notice 
that the periods of these oscillatory motions take values 
in the range between 0.5 and 1.5 ps. This means angu- 
lar frequencies between 4 and 12 ps~^. These numbers 
may be compared with the Rabi frequency, g s/h = 1.5 s 
ps~^, which takes values around 2-4 ps~^, and the char- 
acteristic frequency of Coulomb interactions, fi/h ~ 10 
ps~^. 

Let us stress that we are pumping both single-pair lev- 
els at the same rate, P. The qualitative features of the 
Nstates = 2 problem, mentioned above, do not seem to 
rely on the specific pumping scheme used. The qual- 
itative features dot not change either when relaxation 
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FIG. 7: Number of photons as a function of the photon energy 
for the Nstates ~ 3 dot and cavity parameters (a) P = 3 ps~^, 
>c — 1 ps~^, and (b) P — 1 ps~^, x — 0.1 ps~^. 



from the higher-energy level to the lower-energy one is 
included. 



The iV, 



= 3 dot coupled to aj+i light 



In general, the resonance curve, Nph vs A, in the 
Nstates = 3 system exhibits three peaks, as one can ex- 
pect (Fig. [7^). For certain values of the system parame- 
ters, two or more peaks may overlap, or additional sharp 
peaks (the result of interference) may emerge (see, for 
example. Fig. U}p)- 

Concerning the relative phases between the polariza- 
tion functions and the radiation field, 4> — (f>i, one can 
roughly say the following. In the uppermost resonance in- 
terval, the three single-pair levels cooperate. This means 
that they are all emitting simultaneously, or periodically 
absorbing-emitting almost "in phase" . In the central in- 
terval, in general, there are two levels cooperating, and 
the third is not. For example, two of them are emit- 
ting, and the third is absorbing. Finally, in the lowest 
interval, when one of the levels is emitting, the second 
is absorbing, and the third alternates between the first 
two. Curves are qualitatively similar to the Nstates — 2 
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FIG. 8: (Color online) The Nstates = 2 dot with spin-flip 
processes included: (a) = 5 ps, (b) Ts = 100 ps. See 
explanation in the main text. 



case, and will not be shown. 



D. The Nstatca = 2 dot coupled to different photon 
modes 

In this section, we consider the possibility of spin-flips 
between exciton states with a characteristic rate 2tt/ts. 
Typical values of Ts are around 5 ps or larger—. Then, 
there is a new nontrivial situation for the Nstates = 2 
dot, in which one electron- hole pair is coupled to the 
CT(_|.) photon mode, and the second one to the CT(_) mode. 

The spin-flip term will be added to the equations for 
the occupations, Eq. (|17() . At _B = 7 Teslas, the Zeeman 
splitting between exciton states is around 0.7 meV. The 
lowest exciton state is coupled to the cr(_) photon. Then, 
in the equation for pii (the exciton coupled to the cr(+) 
photon), we add a term 

--Pu(I-Pit), (29) 

and for pi| a term 



— Pu(I-Pit)- (30) 

Notice the Pauli blocking factors in them. Spin-flip tran- 
sitions from pif to pii would require an activation mech- 
anism to overcome the energy barrier, and will not be 
included. 

The results for the system with cavity parameters 
K ~ 0.1 ps~^, P(+) = 5 ps~^, P(_) = 2 ps~^, are shown in 
Fig. IHl The lower figure corresponds to Ts = 100 ps, i.e., 
when the CT(+) and cr(_') dynamics are independent. The 
upper one, to = 5 ps, that is, a rate 27r/Ts w 1.2 ps~^. 
Step-like variations in one of the photon populations are 
transmitted to the other photon number. Interesting 
enough is the region 0.4 < A < 1.3 ps^^, where there are 
only CT(_) photons in spite of the fact that P(+) > P{~)- 

IV. CONCLUDING REMARKS 

We have explicitly written and numerically solved the 
evolution (Bloch-Lamb) equations describing the dynam- 
ics of a quantum dot-microcavity system. The main sim- 
plifications contained in our equations are the following: 
(a) they are mean field equations, i.e., dot not contain 
particle-particle correlations, and (b) the photon field is 
coherent. Assumption (b) is closely related to (a). The 
reward from this simplified approach is that the system of 
equations is relatively small, allowing a qualitative analy- 
sis of its solutions and the extension to larger dots, which 
may support more than one electron-hole pair. 

In the simplest system with Nstates = Ij the qualitative 
analysis revealed the existence of resonance intervals for 
the parameters, where the coupling to the photon field 
of the cavity is optimal and the number of photons in 
the stationary state reaches a maximum. There is also a 
second stable stationary state, which is not reached be- 
cause the used (null) initial conditions is not in its basin 
of attraction. The two stable solutions should be further 
analysed focusing on a possible optical bistability in this 
system.^^. Outside the resonance interval, one of the sta- 
ble states dissapears (suggesting a bifurcation) leading 
to the appearence of stable periodic orbits in which the 
physical magnitudes oscillate. In this oscillatory states, 
the quantum dot is periodically absorbing and emitting 
photons from and to the cavity. Under strong pumping, 
the transition resonance-out of resonance may be almost 
step-like, a fact which may also have interesting applica- 
tions. We have left many questions unexplored such as, 
for example, the width of the resonance interval in A, or 
the transient time before the stationary state is reached, 
as functions of the system parameters. 

Let us stress that the existence of two stable station- 
ary solutions are intuitively related to the two possible 
couplings between an exciton and a photon (the two po- 
lariton branches in a quantum well under weak coupling, 
for example). However, we can not yet understand the 
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way they may arise from the system of hnear equations of 
Ref. |2- We think that the truncation fo the photon Fock 
space is responsible for this difference with our results. 

Larger systems exhibit even more complex dynamics, 
which seems to strongly depend on the position of the 
single-pair levels. Most of the new features are already 
seen in the Nstates = 2 dot, when both pair states are 
coupled to the same photon mode. We showed results 
for a dot in which the energy levels are 8-10 meV ap- 
part, in such a way that individual resonance intervals 
are resolved. Coulomb interactions are shown to have 
drastic effects, not only on the level positions, but also 
causing interferences and altering the very nature of the 
stationary states. Oscillatory states are reached almost 
for any values of the system parameters, even inside the 
resonance intervals. Cooperative or non cooperative be- 
havior of the single-pair levels are obtained in different 
parameter regions. In the Nstates — 3 system, these fea- 



tures are reinforced. 

The account of relaxation between electronic states or 
a change in the pumping scheme used to feed the dot are 
shown to have small impact on this qualitative and semi- 
quantitative picture. Other effects such as the coupling 
to the delocalized states of the wetting layer surrounding 
the self-assembled dot are to be considered. This cou- 
pling to the wetting layer is crucial in other contexts^^. 

Work along some of the abovementioned questions is 
currently in progress. 



Acknowledgments 

The authors acknowledge the Comitee for Research of 
the Universidad de Antioquia for support. 



^ D. Meschede, H. Walther, and G. MuUer, Phys. Rev. Lett. 
54, 551 (1985). 

^ J.I. Perea, D. Porras, and C. Tejedor, 

http: / /arxiv.org/cond-mat /0310570 
^ J.M. Gerard, B. Sermage, B. Gayral, B. Legrand, E. 

Costard, and V. Thierry-Mieg, Phys. Rev. Lett. 81, 1110 

(1998). 

* G.S. Solomon, M. Pelton, and Y. Yamamoto, Phys. Rev. 

Lett. 86, 3903 (2001). 
^ M. Pehon, C. Santori, J. Vuckovic, B. Zhang, G.S. 

Solomon, J. Plant, and Y. Yamamoto, Phys. Rev. Lett. 

89, 233602 (2002). 
® O. Benson, C. Santori, M. Pelton, and Y. Yamamoto, 

Phys. Rev. Lett. 84, 2513 (2000). 
^ T.M. Stace, G.J. Milburn, and C.H.W. Barnes, Phys. Rev. 

B 67, 085317 (2003). 
® J.M. Gerard, D. Barrier, J.Y. Marzin, R. Kuszelewicz, L. 

Manin, E. Costard, V. Thierry-Mieg, and T. Rivera, Appl. 

Phys. Lett. 69, 449 (1996). 
® A. Yariv, Optical Electronics (Saunders College, San Fran- 



cisco, 1991). 

G. Bastard, Wave mechanics applied to semiconductor het- 
ero structures (Les editions de physique, Les Ulis Cedex, 
1998). 

H. Haug and S.W. Koch, Quantum Theory of the Optical 
and Electronic Properties of Semiconductors (World Scien- 
tific, Singapore, 1994). 

M.O. Scully and S. Zubairy, Quantum Optics (Cambridge 
Univ. Press, Cambridge, 2001). 

E. Ott, Chaos in Dynamical Systems (Cambridge Univ. 
Press, Cambridge, 2002). 

J.M. Kikkawa and D.D. Awschalom, Phys. Rev. Lett. 80, 
4313 (1998). 

^5 M. Gurioh, L. Cavigli, G. Khitrova, and H. Gibbs, Phys. 
Stat. Sol. (a) 201, 661 (2004). 

Q.Q. Wang, A. MuUer, P. Bianucci, E. Rossi, Q.K. Xue, 
T. Takagahara, C. Piermarocchi, A. H. MacDonald, and 
C.K. Shih, ,http://arxiv.org/cond-mat/0404465) 



